Directional emission of single photons from small atomic samples 
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We provide a formalism to describe deterministic emission of single photons with tailored spatial 
and temporal profiles from a regular array of multi-level atoms. We assume that a single collective 
excitation is initially shared by all the atoms in a metastable atomic state, and that this state is 
coupled by a classical laser field to an optically excited state which rapidly decays to the ground 
atomic state. Our model accounts for the different field polarization components via re-absorption 
and emission of light by the Zeeman manifold of optically excited states. 
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INTRODUCTION 



Single photons may serve as flying qubits to commu- 
nicate between registers of stationary, material qubits in 
quantum computing architectures and they may be 
applied in protocols for quantum cryptography. In these 
protocols transmission losses over long distances can be 
counteracted by transfer of the light state qubits to quan- 
tum repeaters for purification and entanglement distilla- 
tion [2]. Candidates for stationary qubits that can effec- 
tively interact with single photons are optically thick en- 
sembles of atoms [3-7], rare-earth ions in crystals [sl-Hl], 
vibrational excitations in diamond crystals [12| , as well as 
systems with fewer particles using optical cavities to in- 
crease the interaction with the photon field [l3l-[l7| . The 
systems mentioned can provide a deterministic coupling 
of the material system to a suitably tailored spatial and 
temporal photon wave packet. There are also a number 
of probabilistic protocols, where measurement processes 
herald the successful generation of non-classical excita- 
tions of either the photon field or the medium [l8l-[20| . In 
this article we utilize the fact that an ensemble of atoms 
of just few hundred atoms may interact strongly with a 
single mode of light with a specifically chosen mode func- 
tion. We identify this mode function by calculating the 
emitted field from the atomic ensemble, prepared in a 
collectively excited state. By a time reversal argument, 
the complex conjugate of this emitted field may be in- 
jected on a ground state atomic ensemble and will then 
be fully absorbed at a definite instant of time [2l|, l22j . 
We discuss the possibility to shape the temporal profile 
of the emitted photon, and, in particular, the creation 
of time- symmetric photon wave packets, as such packets 
can then be emitted by one ensemble and absorbed by 
another one in a fully deterministic manner. We focus on 
samples of few hundred atoms, distributed over a few mi- 




FIG. 1: Atoms with three excited substates e~ , e and e + , 
a ground state g and a long lived state / are arranged in a 
regular lattice. The state / is coupled to the state e +1 by a 
classical laser field. The excited states decay to the ground 
state emitting a photon with a - , ty and cr + polarized photons, 
respectively. The direction of the emission is determined by 
the wave number of the atomic g — f coherence and the wave 
number of the out-coupling field, k e m — k g f — kh- 



crometer spatial extent. In such ensembles, the Rydberg 
blockade interaction may be used to establish singly ex- 
cited states and, subsequently, single photon states [23l — 
[26[, while photonic qubits, collectively absorbed by the 
atoms may be manipulated by Rydberg state mediated 
quantum gate operations [2l|, [27|, [28[ . 

The collective interaction of light with ensembles of 
absorbers and scatterers has been an active field of study 
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since the early days of electromagnetism, while collective 
phenomena in spontaneous emission received wide atten- 
tion with the pioneering work on Dicke superradiance 
from population inverted samples [29]. Early studies of 
collective emission from ensembles with few excitations 
(30l-[34j (see also [35] and references therein) have been 
followed by a recent flourishing of analyses [2l|, l36l - [43j , 
which apply a Born-Markov approximation and eliminate 
the field degrees of freedom to obtain coupled equations 
for the atomic excited state amplitudes. Approximate 
solutions to these equations may be derived, e.g., with 
the assumption of a scalar description of the field, but 
for only few hundred atoms, they may also be solved di- 
rectly on a computer. 

In this manuscript we generalize the previous analyses 
to account for the full vector character of the quantized 
radiation field. We establish coupled equations for ex- 
cited state amplitudes on a suitable set of atomic Zee- 
man sub- levels, emitting and reabsorbing the different 
polarization components of the field, and we solve the 
equations numerically to identify the full temporal, spa- 
tial and polarization content of the emitted light. 

In Sec. II we derive the coupled atomic equations un- 
der the Born-Markov approximation. In Sec. Ill, we 
present numerical results for the photon modes emitted 
by samples of atoms with different spatial geometries. In 
Sec. IV, we describe the use of a coupling laser field to 
control the temporal shape of the emitted photon wave 
packet, and in Sec. V, we present a brief conclusion and 
outlook of the work. 



II. DIPOLE-DIPOLE INTERACTION 

We want to describe the experimental situation where 
a collection of N atoms can be prepared in a single ground 
state g, and where a suitable, symmetric excitation mech- 
anism allows the preparation of a state 



N 
3 = 1 



a>j\gig2-,fj,-gN) 



(i) 



with a single atom transferred to the metastable state 
/, see Fig. 1. We assume that N atoms are located at 
the positions fj (j = 1,..^). With plane wave excita- 
tion laser fields, the amplitudes ctj may have have equal 
magnitude and they depend on the phase of the fields 
at the atomic locations, aj = ^= exp(— ik g f • fj). The 
Rydberg blockade mechanism may restrict the system to 
precisely one excitation, and the state /, sketched in the 
figure, may indeed represent a long-lived Rydberg state, 
or a Raman process via a Rydberg state, may transfer 
precisely one atom to a long-lived low lying atomic state. 

To release a photon from the system, we use a classical 
laser field with the Rabi frequency VLl to drive the atomic 
/-state amplitude into an optically excited state e, with a 
strong dipole coupling to the ground state g. The system 
now acts as an antenna array for dipole radiation on the 



transition e — g, and this is the cause of the desired direc- 
tionality of the emitted light. As indicated in Fig. [TJ the 
initially populated states may be extremal Zeeman sub- 
levels with well defined polarization selection rules, and 
the photon emitted on the e — g transition may be <j + 
polarized with respect to the atomic quantization axis. 
This field, however, may be reabsorbed by another atom 
located in an arbitrary direction from the emitter, and 
here, the expansion of the field on polarization compo- 
nents permits excitation with selection rules Am = 0, ±1. 
To describe the many-atom emission, we thus have to 
include other Zeeman sublevels than the ones initially 
populated. This motivates the model depicted in Fig. 
1, with unique states g and /, and three excited states 
e°, e ±x , corresponding to a J = — J = 1 optical tran- 
sition. This configuration is the simplest extension of a 
two-level model atoms that allows us to fully take into 
account the polarization of the emitted and re-absorbed 
light as well as the resulting dipole-dipole interactions 
between the atoms. 

In the following we will use the short hand notation 
for singly excited states of the atomic ensemble, \fj) = 
|#i#2---/j---<7iv)> and similarly for \e v ), with v = 0, ±1. 

In the dipole approximation the interaction of atoms 
with photons is described by a Hamiltonian (3lj : 



H = Hq + Hi ri 
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where 
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is the atom-field Hamiltonian and the interaction part is 

Hint = Hl + Hy • (4) 

The semi-classical coupling to the initial long-lived state 
is 
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where el is the polarization direction of the coupling field 
with the optical frequency ujl and <xj> e = df e \fj)(e~j~ 1 \. 
The direction of the dipole moment for this transition 
df e we further assume to be parallel to e^, so that the 
transfer of amplitude happens exclusively to the state 
le+i). 

The coupling of the atomic dipole between \g) and \e) 
to the quantized radiation field modes is described by 
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Here, the atomic dipole operator is defined as a J ge = 

^2t = _id J gu \g)(e 1 j\, where dP gv is the unit vector in the 
direction of the corresponding dipole moment for the 
\g) — \e v ) transition. a^ x is the annihilation operator 

of a vacuum electro-magnetic field mode k with the po- 
larization A in the direction e^ A . The atom-field coupling 

strength is gk = d eg ( 2 ^y k ) X ^ with the quantization vol- 
ume V and uok = cfc. 

We henceforth ignore spontaneous emission on the e- 
/ transition; this may on the one hand be chosen as a 
transition with a weaker dipole moment, and on the other 
hand it does not experience the collective enhancement, 
that we shall observe on the e-g transition. 

We expand the time dependent solution of the 
Schrodinger equation for N atoms and the field as a su- 
perposition of Fock states with a single atomic or pho- 
tonic excitation 

N 
3=1 

+ E E ^We-^ t |0)|ej}|0)+ 

j=l v=-\ 

+EE^We- iWfct io)i5)ii fe - A ), (7) 

h A 

where \g) represents the state with all atoms in the 
ground state. 

Note that we use the rotating wave approximation 
(RWA) in Eq. (j5]) because we treat this transition semi- 
classically. In the quantized atom-light interaction, 
Eq. ([6]) , the RWA shows a rather intricate interplay with 
the role of virtual photon processes and a seemingly coin- 
cidental equivalence of terms as discussed for the case of 
two- level atoms in [3l|, (36[ . This holds for the terms de- 
scribing atom-atom interactions. A generalization of this 
non-trivial discussion to the case of multilevel atoms in- 
teracting with a quantized vector field will be presented 
elsewhere Q. This equivalence of terms allows us to 
use the RWA in our treatment of quantized atom-light 
interaction in Eq. (|6j). 

Substitution of Eq. ([7]) with the initial condition 
e^ A (0) = into the Schrodinger equation with the Hamil- 
tonian from Eq. ([2]) yields the formal solution 

N +1 
3=1 v=-\ 

x J dr/3J(r)e- i(u;o -^ )r . (8) 

Using the Markovian approximation [30| the atomic co- 
efficient /3j (r) can be approximated by /3j (t) in Eq. J5J) 
and taken outside the integral. This allows us to substi- 
tute the photon amplitudes by expressions involving only 
atomic state amplitudes, which thus obey a closed set of 
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FIG. 2: A snapshot of the total intensity emitted by a cloud 
of 14 x 14 x 10 (blue) and 3 x 3 x 10 (green) atoms with the 
lattice spacing d = 0.25 Ao and a single atom using a square 
out-coupling pulse (dotted line). Due to the coupling between 
different decay channels in the system, we observe beat-like 
behavior of the emitted light. See as well Sec. IIHI for more 
details. Parameters used for this simulation are Ql — 8.2 T, 
the single atom detuning of Ql light is and the pulse length 

t w = 0.2 r _1 . 



equations: 

hi = ^c^/--^)*ft +1 , (9) 




- \ E E ( J - (4 • ^ • 4*) ( 10 ) 

j=l v=-l 

Here V is the single atom decay rate from |e) to and 
A Lamb is the single atom Lamb shift. This term contains 
an infinite integral, where a suitable cut-off should be 
applied to yield a finite physical value [45|. After this 
procedure this shift can in principle be absorbed into 
the definition (the measured value) of the energy of the 
atomic state \e) in Eq. (0. 

The second rank tensor 
i*<?(koRij) with 

fi kR) = -(l-RR)-L r l + 

+ 2l J - mR ) \~Wf-^Rf)' (H) 
_ N 3 AA \ cost kR) 

--[I -3RR) [-J^T + ^r) ( 12 ) 
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and Ri j = fl — fj accounts for the field mediated inter- 

action between the atoms / and j. Here / is the unity 
tensor and RR is the projection onto the direction given 
by R H- 

Eqs. ([9])- ([T0]) account for the zero-photon subspace 
component of the total wavefunction of the atoms and 
the quantized field, 

N N +1 

l*o(«)> =$>(«)!/.,•> + £ £ ^(t)\ej), (13) 

j=l j = l v=-l 

and this state component is described by an effective non- 
Hermitian Hamiltonian 

with 

jfjeff = _ ih _ iALamb \ J2 £ |e?>(e?| (15) 
vz J i=i ,=-i 

describing individual single atom effects, 
so N 

^ = -^E^-^l/^Xeri (16) 

Z = l 

giving the coupling from the long-lived state / and 

H e JJ = in T -Y. E (i-Mx 

j,l=l u,rj=—l 

x(4-^-dj v )|e?><e?| (17) 

expressing collective dispersive and dissipative effects be- 
tween the atoms in the ensemble. 

In the rest of this paper we are interested in the spa- 
tial and temporal emission profiles. The V-part °f ^ 
defines the Hermitian part of H^JJ , i.e. the coherent ex- 
change of excitation between the atoms. The / part 
conversely defines the anti-Hermitian part correspond- 
ing to the decay by emission of light, and thus popu- 
lation of the one-photon quantum state component. If 
the Hermitian part is diagonalized, delocalized orthogo- 
nal eigenmodes with collectively " Lamb shifted" energies 
are obtained. If the anti-Hermitian part is diagonalized, 
delocalized orthogonal independently decaying modes are 
obtained. Some of these modes have a decay time longer 
than T, i.e., they are Dicke subradiant modes [29], and 
some decay faster than T, i.e., they are Dicke superra- 
diant modes, see [2l| . (36j and references therein. The 
Hermitian and the anti-Hermitian parts do not commute, 
and hence the eigenmodes of the full Hamiltonian H e ff, 
which provide the time dependent atomic state as a sin- 
gle sum of complex exponentially weighted vectors, are 
not orthogonal. This will result in coupling between dif- 
ferent decay channels and a quantum beat-like behavior 



for the decay modes [Io[. The blue line in Fig. [2] shows 
a snapshot of the total intensity emitted by a cloud of 
14 x 14 x 10 atoms excited from the levels \fj) to the lev- 
els le^ 1 ) by a square pulse (dotted line). The initial 
emission rate is faster than the single atom emission rate 
T presented by the red line. Additionally, we see a mod- 
ulation of the emitted intensity due to the coupling of 
the different non-orthogonal decay modes. With the re- 
duction of the cloud size, the coupling between the decay 
modes becomes weaker, but the superradiance behavior 
remains pronounced as shown by the green line in Fig. [2] 
for the array of only 3 x 3 x 10 atoms. 

III. SPATIAL PHOTON MODES 

The eigenmode expansion of Eqs. ([9])- ([T0]) for the 
atomic excitation amplitudes formally yields a solution 
for each individual atomic excited state amplitude as a 
sum of exponential functions of time with complex ar- 
guments. Once the atomic evolution is determined, the 
light emission is given by the integrals in Eq. (j5J). We 
note that for the relevant time scales, longer than the en- 
semble excited state lifetime, these integrals involve only 
decaying exponential functions, and they may in prac- 
tice be extended to infinity. This allows to determine the 
(far) field eigenmode expansion coefficients as algebraic 
expressions involving the mode expansion coefficients di- 
vided by the sum of the complex eigenvalues and the 
frequency difference appearing explicitly in Eq. (|8j) . This 
means that we can readily determine the field amplitudes 
on any chosen set of field modes after diagonalization of 
the atomic problem, at a cost that depends only on the 
number of atoms, and calculations with even thousands 
of atoms are realistic. 

The probability to detect a photon at a position r at 
time instance t, much later than L/c, where L is the 
linear sample length, is given by |46| 

h (r,t) = (m\4-H^4 +) (r)\m). (18) 

Here e is the handedness of the photon, i.e., its circu- 
lar polarization along the line connecting the atomic en- 
semble and the detector at position f. The positive fre- 
quency component of the desired polarization is given by 

£e (+) (r) = (e- 2?(+)(r)) with |l| 

£(+) (F) =iJ2J2 Kg&^fae** ( 19 ) 

q cr 

See Appendix [A] for the full derivation of I e . 

In this section we assume that the atomic system is 
initially prepared in a so-called timed Dicke state 

a,(0) = -^e-^^,/?J(0) = 0. (20) 

This state is coupled to |e +1 ) by switching on the laser 
field described by the Rabi frequency VLl, and we first 
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FIG. 3: Emission of a single photon from a small atomic sample of 3 x 3 x 8 atoms arranged in a lattice with the spacing 
d = 0.60 Ao- a) Integrated intensity as function of the distance from the atomic sample (solid line) and the out-coupling pulse 
shape (dotted line), b) Spatial distribution of intensity of left-handed polarized light for a given time t — 200 Xo/c. The 
intensity scale is given in arbitrary units. The spatial coordinates are given in units of Ao- c) Same as b), but for the opposite 
handedness of the light, d-e) The angular intensity distributions of light at R = 199 c/Y corresponding to Figs, b) and c) at 
the moment of the maximum intensity. Note that the intensity scales for the two polarizations differ by an order of magnitude. 
Parameters used for this simulation are Ql = 2 F and the single atom detuning of the Ql field is 10 T. 



study the case where this coupling field is kept constant. 
In order to avoid resonant coupling to individual modes 
of Eq. ([TT]) , we assume as well that the Q l field detuning 
from the single atom resonance |/) — |e +1 ) is larger than 
the Rabi frequency f^, see Sec. II VI for further discussion. 

We first consider a sample of 3 x 3 x 8 atoms arranged 
in a lattice along the a>, y- and z- directions, respectively, 
with the lattice period d > Xq/2. This case would corre- 
spond to trapping of alkali atoms on the ujo transition in 
a red detuned lattice. 

The calculated emission of a single photon from this 
atomic sample is presented in Fig. [3l Figure [3^ shows 
the intensity of the field, integrated over directions, at 
different distances from the atomic sample. The field 
propagates at the speed of light, and this snapshot of 
the intensity distribution with distance reflects how the 
state \f) population has gradually decayed via the opti- 



cally excited states since the coupling field was switched 
on. The most prominent feature is the overall exponen- 
tially decaying shape, but the figure also show the resid- 
ual beat-like behavior. This beat-like behavior becomes 
even more pronounced with the increase of the array size 
as was shown in Fig. [2] and is in agreement with our 
discussion of the real and imaginary eigenvalues of the 
problem. Figures [3}} and c show the spatial intensity dis- 
tribution, indicated by colors, for the same parameters as 
in part a of the figure. Figure ^jp shows the intensity dis- 
tribution for the polarization of the emitted light, which 
is expected to be dominant for the level scheme by the 
dipole selection rule. We observe that most of the light 
is emitted in a narrow forward peak. Figure [3t shows 
that in the backward direction, a small component (note 
the different scale) is emitted with the opposite polariza- 
tion, again in agreement with the atomic dipole selection 
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rule. To quantify the angular distribution of the emit- 
ted light, Figs. [3]i,e show polar plots of the intensity of 
the left-handed and right-handed light components at the 
distance R = 199 c/T. 

The directionality of the emitted light as a function of 
the lattice spacing was analytically studied by Porras and 
Cirac [39] for the case of two- level atoms. In this model 
the radiative pattern can be factor ized into two parts: 
1) collective scalar field radiative effects and 2) a dipole 
radiative pattern of individual atoms. In a general case 
of multiple exited states and a vector electro-magnetic 
field, atomic reabsorption of photons with a resulting re- 
distribution of angular momentum associated with the 
directions of propagation between the atoms, does not 
justify such a simple factorization of the radiation pat- 
tern and a scalar field approximation. Our calculations 
show, that these reabsorption effects become important 
for the case of strongly coupled atoms, i.e., for atomic 
spacings d < 0.25 Ao, and in the case of frustrated for- 
ward emission even for spacings d « Ao, see Fig. HJ In 
both cases the subradiant modes play an important role. 

Although the detailed analysis of these cases is beyond 
the scope of the current work, we give short example 
for the case of frustrated emission. We consider emis- 
sion by an array of 8 x 8 x 8 atoms with the spacing 
d = 0.60 Ao- The atoms are initially prepared in the su- 
perposition timed Dicke state of the |/) levels Eq. ([20]) as 
in all previous examples, favoring the emission in the z- 
direction, but the storage level |/) is coupled now to the 
short lived state |e°). This state can directly decay to 
the ground state with the emission of a 7r-polarized pho- 
ton, relative to the z-axis, see Fig. HJ Since the angular 
emission pattern of this transition has zero intensity in 
the z-direction, which is the preferred emission direction 
of this atomic sample, we arrive at the conflict between 
the preferred and allowed emissions. An example of the 
frustrated emission is presented in Fig. 2J Here we plot a 
snapshot of the angular distribution at R = 208 c/T and 
the integrated intensity emitted by an array of 8 x 8 x 8 
atoms with d = 0.60 Ao and a long coupling pulse VLl 
(dotted line). Figures H^i-b show the case, where the cou- 
pling of \g) to all |e _1 ), |e°) and |e +1 ) states is allowed. 
Fig. [4]>d show the situation of emission from two level 
atoms, where only \g) and |e°) are coupled. Both the 
angular dependence and the time evolution show differ- 
ences in these two situations. This difference stems from 
the reabsorption of virtual photons with different polar- 
izations. 

Our calculations also show, that the analysis of [39| 
is qualitatively correct for atoms with larger separation. 
The directionality of emission in this case depends on 
the lattice spacing and can be understood as interference 
of Bragg scattering contributions. For a critical spacing 
of d = Ao/2 or an integer multiple of Ao/2, the photon 
is emitted in the forward and backward directions with 
equal probability, see Fig. [5ji and e. With an increased 
value of d this symmetry is broken and the forward 
emission peak becomes dominant, see Fig. [5]>d and f. 



Due to the diffraction-like effects the directionality of 
the emitted light improves with the increase of the array 
size. Figures [5b, g and h show the angular distribution 
of emitted light for lattices with the same spacing, but 
with an increasing number of atoms. 

In this section we have studied the collective emission 
of single photons and we have analyzed details of their 
directional distribution and their associated polarization 
properties. In the next section we turn to the temporal 
shape of the emitted light pulses and to the experimental 
means to control the light mode of the emitted photon. 

IV. TEMPORAL PHOTON MODES 

The angular emission patterns shown in the previous 
section were all calculated assuming a constant intensity 
of the laser which drives the f <-> e transition. This leads 
to a highly asymmetric temporal profile of the emitted 
light shown in Fig. [3^. By controlling the temporal pro- 
file of the coupling fiz,(£) one can control the temporal 
shape of the emitted light, and this can be used to trans- 
fer the atomic excitation from the state |/) to, e.g., a 
temporary symmetric emitted photon wave packet. Such 
a wave packet can be reabsorbed in a second atomic en- 
semble, if one employs the time reversed control field. 
Therefore we may imagine a collection of atomic ensem- 
bles as quantum repeater stations, where photon pulses 
are absorbed and reemitted, possibly after suitable en- 
tanglement distillation and state purification 0, [ll| . 

To design appropriate outcoupling laser fields, we start 
with a system prepared in the states \fj) described by 
Eq. ([20]) with a vanishing coupling Ql. By gradually 
increasing the coupling strength we transfer population 
to the states Je^ 1 ), which decay to the ground state by 
emission of light. The intensity of the emitted light is 
given by the population of the excited state Je^ 1 ), and 
our goal is thus to control this population. 

Due to the complexity of the non-orthogonal eigen- 
modes of the coupling Hamiltonian, see Sec. [Til we fo- 
cus here on the conceptually easiest strategy, which is 
an adiabatic out-coupling. By introducing a detuning 
5 — ujf e — ujl and fi\ = fi\e~ lbt we can rewrite Eq. (|9j)- 

m 

* = ^Pt\ (21) 

k 1 = ^ra l + i(S + A Lamb ) - 
r N 

~2 £ (^-y^-dU)^ 1 - ( 22 ) 

Since the state \fj) is coupled by a strong field VLl to 
le^ 1 ), whereas the substates {ej 1 } and \e®) are only cou- 
pled by virtual photon processes to the state Je^ 1 ), we 
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FIG. 4: Frustrated emission from a sample of 8 x 8 x 8 with the lattice spacing d = 0.60 Ao- a) A snapshot of the angular 
distribution at R = 208 c/T and b) angularly integrated intensity of the emitted light for the case of all |e _1 ), |e°) and |e +1 ) 
levels are included. c)-d) The corresponding graphs, if only |e°) level is included. The coupling laser parameters are the same 
as in Fig. [3] 



have neglected the populations of the {ej 1 ) and \e®) sub- 
states. Assuming further that the detuning S is much 
larger than the Rabi frequency f}/,, the single atom de- 
cay rate T and the Lamb shifts, we adiabatically elimi- 
nate the state /3j +1 

r N 

-*28 £ (^V^-i.)^' ( 2 3) 

We can see directly that the leading term of is of 
the order O (pp), since a\ is of the order of one. There- 
fore the population of the state is always much 



smaller than unity and quickly decays to the ground 
state. Hence the intensity of the emitted light is ap- 
proximately given by the population of the state Je^ 1 ). 
Substituting Eq. ([23]) into Eq. ([2T]) and neglecting terms 

of order higher than O (jfir^ we arrive at a closed set of 

equations for the a coefficients: 

ai = -i Ang ht cli - ~l^ a i~ 

T ! N 

-y E ( 24 ) 

with the light shift Au g ht — and the effective decay 
rate V = Tj^ of the metastable state \fi) . 
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FIG. 5: Total angular distribution of emitted light for a lattice 3x3x8 atoms for the cases a) d = 0.50 Ao, b) d = 0.53 Ao, 
c) d = 0.60 Ao, d) d = 0.90 A , e) d = 1.00 A , f) d = 1.10 A . The total angular distribution for lattices g) 6 x 6 x 8 and h) 
8x8x8 with d = 0.60 Ao. The coupling laser parameters are the same as in Fig. [3] 



Similar to Eq. ([22]) , the first line of this equation de- 
scribes the dynamics of an isolated atom. The second 
line, which depends on the geometry, is the contribution 
from the interactions with all other atoms of the sample 
via virtual photon exchange processes. Since the collec- 
tive contributions both in Eq. ([22]) and in Eq. ([M]) have 
the same form, and the temporal evolution of the popu- 
lation follows ai as 



5+1 

# = 2S ai > 



(25) 



up to the order O (^p), we conclude that the spatiotem- 
poral mode function occupied by emitted photon has the 
same structure, up to a radial scaling factor, as in the 
case studied in Sec. IIIII 

In order to tailor the temporal dynamics of the emitted 
light, we now allow temporal modulation of the control 
field Q L (t) = Q L of(t) with < f(t) < 1 for all times. 
With this notation Eq. ([24]) becomes 



1 da,(t) _ A Tn lo n(f) 
JW dt ~ ai[t) ~ 2# ai{ty 



~\% £ (26) 
In contrast to a resonant outcoupling case described by 



N 



Vasilev et.al. [49], in our adiabatic outcoupling case there 
is no general analytic solution of this equation. Never- 
theless we can find a connection between the solution 
afi(i) of Eq. ([26]) with a constant Ql, i.e. / = 1, and the 
solution ai(t) for an arbitrary /(£). Defining 



<t)= ff(t') 2 dt', 

JO 



we observe that 

obeys the equation 
dh(T) 



bi(r(t)) =ai(t), 



(27) 



(28) 



dr 



4 *w 

N 



r o 2 

2 4^ bl{T) 



3 = h&l 



Since this equation is identical to Eq. ([26]) for constant /, 
we directly obtain its formal solution bi(r) = a°(r), and 
hence, the general solution to Eq. ([26]) reads 

(30) 



a l {t) = aP l {r{t)). 



I.e., the time evolution of ai(t) can be advanced or re- 
tarded with respect to a® (t) by choosing the appropriate 
function f(t) in ([ZZD- 
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a) b) c) 




100 200 300 t x(t) 100 200 300 

R < c/r ) time (a.u.) R (c/r) 

FIG. 6: Temporal shaping of the photon emission from an atomic array of 3 x 3 x 8 atoms with d = 0.60 Ao- a) Output intensity 
(blue dashed line) using a constant pulse (thin, red dashed line) with Qlo = 10.5 Y and 5 = 120 T. With a specially designed 
outcoupling pulse with peak Rabi frequency Qlo = 42.0 Y and the temporal shape f(t) (thin red line), we obtain a Gaussian 
temporal shape at R = 260 c/Y of the emitted light intensity (blue line), b) Integrated intensities n°(t) (dashed) and n(t) (thick 
solid curve) and a graphical solution (red arrows) of the equation n(t) — n°(r(t)). c) Outcoupling of a double peaked pulse 
(thick blue line) using Q^o = 38.85 Y and f(t) represented by the thin red line. 



To calculate the temporal pulse shape f(t) leading to 
any desired output intensity I(t), we observe that the 
number of photons emitted under constant amplitude 
driving is related to the population in the initial atomic 
state 

nt N 

n°(t)= / dt'I°(t')) = l-J2\*Ht)\ 2 : (31) 

J ° 1=1 

while the corresponding number of photons in the desired 
field is 

n(t)= / <tf/(0 = l-y>(*)| 2 - (32) 
Jo i=i 

Eq. ([30|) implies the equation for r(t): 

n(t)=n°(r(t)). (33) 

Fig. [6^ shows the emitted intensities (dashed blue curve) 
obtained with a constant outcoupling amplitude (red, 
dashed curve) and a Gaussian pulse with 1/e width 
15 T _1 at R = 260 c/Y (dark blue curve) obtained 
with the calculated amplitude f(t) (thin red curve), re- 
spectively. The corresponding integrated intensities are 
shown in Fig. [6b, which also illustrates the numerical 
procedure to solve (|33|) and establish the correspondence 
between the values r(t) and t, as indicated by the red 
arrows. By numerically differentiating the function r(t) 
and using Eq. (|2T|) we directly obtain the temporal pro- 
file of the outcoupling pulse f(t) leading to the desired 
I if). This function is shown as a red curve in Fig. [6^i; the 
spikes are due to precision errors in the calculation of the 
derivative of r(t). We insert this outcoupling pulse into 
our simulation of the full coupled equations and we calcu- 
late the resulting emitted intensity which is, indeed, the 
result, presented as the blue line in Fig. [6^. To demon- 
strate our ability to couple out arbitrary pulse shapes, 



we aim in Fig. [6b at a double peaked pulse with equal 
heights at R = 160 c/Y and R = 260 c/Y. The corre- 
sponding profile of f(t) is shown as a red line in Fig. [6b, 
and again, the solid, blue curve shows the outcome of the 
simulation of the full set of coupled atomic equations. 

We presented here a recipe for photon shaping starting 
from a constant amplitude control pulse. Using this as a 
reference, it is possible to relate the outcome of different 
control pulses by a suitable parametrization of the time 
argument, and to use the above method with experimen- 
tally measured intensities. While we assumed a varying 
amplitude but a constant phase of the control field f2i,(t), 
to produce more complicated single photon wave packets, 
it may be worthwhile to study also complex valued /(£). 



V. CONCLUSION 



We studied correlated spontaneous emission from small 
arrays of atoms. Our atomic model is the simplest gen- 
eralization of a two-level atom, which allows to fully ac- 
count for polarization effects. Numerically solving the 
equations of motion we demonstrated the possibility to 
deterministically generate single photons with a well de- 
fined spatial emission profile, which can be controlled 
by the geometry of the array. Finally, we presented a 
simple method for temporal shaping of the emitted pho- 
ton, which is a prerequisite for high fidelity interfacing 
photon wave packets to atomic systems. Specifically we 
presented a recipe for generating symmetric photon wave 
packets, as they can be efficiently re- absorbed by another 
similar atomic sample. 
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Appendix A 

We derive the explicit dependence of I e (r, t) on the 
atomic part 3 of the wave function. We rewrite first 
Eq. (HHD as [46|: 



J e = ^|^-)|0)b)|0)(0|^|(0|^ + )|^). (Al) 
The photon "wave function" is 



|* e ) = <0|<fl|(0|25i + > (f)W*)> = 



i(k-r-u k t) 



(A2) 



with e^ A from Eq. ©. 

As the second step, we simplify this expression. The 
summation over polarizations is performed using a recipe 
from the Appendix of Smith et.al. [48|: 



N 1 



k 3=1 ^=-1 



X %fe e 



2 i(k-Rj-uj k t) 



f dT^(T) t 

Jo 



-i(ujQ—uj k )i 



(A3) 



with Rj = f — fj. After replacing summation over the 
wave- vectors J2k by an integral ^ 2 ^ 3 J °° duJkU 2 J dQ(k) 
we get 



x / dr^(r) x e -*(<"o-<-*)T ; 



x y dn(jfe) (e • (V - jfejfe) • e iE ^ . (A4) 

The angular integration can already be performed at this 
step |42j resulting in 



N 1 



1=1 l/=-l V ' • y ° 



x / cZr/3^(T)e- i(w °- Wfc)T x 



4tt (V- ?(^) • 4,) (A5) 



with C 



Therefore the photon intensity up to the order of 
O ((kR)~ 2 ) with r >> rj is 



j = l L/=-l ,/U 

x / &/c 2 [e 



dTpUT)e- ickoT x 



ikRj-\-ick{r — t) 



ikRj+ick(r—t) 



(A6) 



where 



have 



introduced 



s 8 ^2^ c p 9 ■ (V - JSjiij) • . Since during the 
emission the value of k is peaked around the atomic 
resonance ko where the last time integral is relevant, we 
can replace k 2 with k% and extend the lower integral limit 
to — oo [ig, i.e. the Weisskopf-Wi gner approximation. 
Using the definition of i5-function 5(t) = ^ f^° oQ dke lkt , 
Eq. (|A6|I becomes 



N 



x [J dT^(T)e- ickoT 5 | i 



(A7) 



Since the last integral with the delta function is always 
zero, we arrive at 

j=i v=-i 3 ^ c ' 

x (e • (Y - RjRj) ■ dQ . (A8) 

Finally, we arrive at the desired explicit dependence of 
the photon intensity on the atomic part f3: 



/,(f,t) = (*«!*«> = ^x 



ikoff-irj-rj/)) , 



X 



(*■(?■ 



(A9) 



where we have used « r for /3 and in the denominator, 
but kept the significant term Rj = r — (f • fj ) in the 
exponents. 

This equation allows a very transparent physical inter- 
pretation for the case of many noninter acting atoms. We 
suppose the states \fj) are fully mapped onto the cor- 



responding le^ 1 ), i.e. an- 
atomic evolution is described by 

1 



in Eq. (O. Therefore the 
? 'e"5 t , (A10) 



//.-, 
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$(*) = o, 

and 

which is the solution of Eq. ([TO]) for koRij ^> 1, i.e. 

vanishing coupling ^ij. The corresponding emitted in- 
tensity is then 

J = l 

^2 l4 AT TV 

j=l j'=l 

(All) 

with 



= (e • (V - 4-4) • i^) (e • (V - R r R r ) ■ d^) . 

(A12) 

Equation (|A12|) can be further simplified using = 
ff(l + 0(^f)), and assuming that all dipoles are polar- 
ized in the same direction, i.e. d J 9i +i = I n this 
case the function Cjj has an interpretation as the an- 
gular dependence of the dipole emission pattern for a 



given helicity e expressed in a tensor form. In partic- 
ular, if we introduce spherical coordinate system along 
the dg^o direction and sum over all polarizations, we ar- 
rive at the well known dipole pattern angular dependence 
^(l-\-cos 2 0) for a dipole emitting on the |e +1 ) — > \g) tran- 
sition. 

Correspondingly, in the direction r = we have 

/.Wtl^x^e-^, (A13) 

while for all other directions r ^ the exponents in 
Eq. (|A11|) average out and give 

M^)«lxfe- r H)q. (A14) 

i.e. 1/N 2 reduced emission intensity. This is the well 
known mechanism of directed emission from an atomic 
sample. Note as well that in this example, and in the 
more general case of spin polarized two level atoms, i.e. 
no coupling of v — +1 states to the other v = 0,-1 
states, the single atom dipole emission pattern factorises 
out in the intensity formula Eq. (|A9|) as was pointed out 
by Porras and Cirac [39| . Nevertheless this is not true in 
the general case considered here due to coupling of the 
v = 0,-1 levels to the level v — +1 via virtual photons 
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